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ABSTRACT 

A new analysis of velocity distribution, optical photometry and X-ray surface 
brightness from ROSAT PSPC data of the galaxy cluster A 2319 is presented. The 
temperature profile derived from ASCA data ( [Markevitch 1996] ) is taken into account. 
A method to check the hydrostatic model in the presence of a temperature gradient 
is proposed. Consistency of the hydrostatic isothermal model and the explanation 
of the "/3-discrepancy" are discussed. Galaxy and gas density profiles of the main 
component A 2319 A are derived, allowing for the effect of the secondary component 
A 2319 B. The inadequacy of a polytropic model, which would produce a binding mass 
decrease with respect to the isothermal /?-model, is discussed. A simple interpolation 
of the temperature profile provides instead an increase of the binding mass and a lower 
baryon fraction thus mitigating the "baryon catastrophe". Assuming as typical the 
value ft, ~ 0.2, a comparison with the most recent estimate of Q^ ucl implies for the 
cosmological parameter Q a ^ 0.4. 

Subject headings: dark matter - galaxies: clusters: individual (A 2319) - galaxies: 
photometry - intergalactic medium - X-ray: galaxies 



1. Introduction 

Detailed studies of the matter distribution in clusters of galaxies provide important clues on 
the growth of condensations and the evolution of the Universe. From X-ray observations it is 



possible to derive both the gas and the total binding mass distributions, under the assumption of 
hydrostatic equilibrium. Optical data, i.e. galaxy photometry and redshifts, combined with X-ray 
observations allow to check the validity of the hydrostatic equilibrium assumption and derive the 



spatial distribution of the dark matter. Most analyses in the past (pones and Forman 1984 




Cowie, 


Henriksen and Mushotzky 1987 


. Hughes, Gorenstein and Fabricant 1985, Hughes 1989], Gerbal 





et al. 1992, Briel, Henry and Bohringer 1992 , Durret et al. 1994| , David, Jones and Forman 
1995, Cirimele, Nesci and Trevese 1997| ) were based on the further simplifying assumption that 



the gas is isothermal, at least within about 1 hg Mpc, possibly with the exclusion of a central 
cooling flow region (see e.g. Fabian, Nulsen and Canizares 1984, White, Jones and Forman 1997). 
This leads to the /3-model ( Cavaliere and Fusco-Femiano 1976| ) which predicts that the dynamical 



parameter f3 spec = ^™pf r , representing the ratio between the energy per unit mass of galaxies and 
gas, equals the morphological parameter (5 fa defined by the fit of the gas density distribution with 
a King profile. The observations show that on aversige ft fit ^ ftspec flS&rctziD 1986| , |Evrard 1990| ). 
However Bahcall and Lubin (1994) ascribed this " /3-discrepancy" to the underestimate of the 
slope of the galaxy density profile, appearing in the hydrostatic equilibrium equation, rather than 
to a failure of the model. In their X-ray-Optical analysis of a sample of Abell clusters, Cirimele, 
Nesci, & Trevese (1997) (CNT) found that logp gas = j3xo ^ogPgai + C in a wide range of densities, 
as predicted by the hydrostatic isothermal equilibrium. This allows to define, for each cluster, 
a morphological parameter flxo-i independent of any analytical representation of p gas and p ga i- 
A comparison of (3xo with f3 spe c supports the explanation of the " /3-discrepancy" suggested by 

at least for several galaxy clusters 



Bahcall and Lubin 1994| and the consistency of the "/3-model'' 
of regular and relaxed appearance. The gas and binding mass distributions thus derived provide a 
typical value of the baryon fraction fs, of the order of 0.2 within 1-2 h^ 1 Mpc ( pirimele, Nesci and 
[Trevese 199% |Evrard 19971 , [Ettori fc Fabian 1999| ,Mohr et al., 1999). This relatively high value, 
compared with the results of primordial nucleosynthesis calculation ( |Walker et al. 1991 , Olive 



[Steigman 1995 ,but see |Burles and Tytler 1998 ) implies that either the cosmological parameter 
is smaller than unity ( White et al. 19~93|) , or fs is not representative of the cosmic value, and 



galaxy clusters are surrounded by extended halos of non baryonic dark matter (White and Fabian 



1995). The latter hypothesis raises the problem of understanding the mechanisms of a large scale 



baryonic segregation. 

However, recently ASCA data have provided direct evidences of gas temperature gradients in 



the outer regions of several 


galaxy clusters (Arnaud 1994, 


Markevitch et al. 1994, 


Markevitch et 


al. 1996, 


[kebe et al. 1996, 


Markevitch 1996|, 


Markevitch et al. 1998 


) . According 


to Markevitch 



1996 , in the outer regions some clusters show a polytropic index even greater than 5/3, which is 



inconsistent with the hydrostatic equilibrium conditions. According to Ettori fe Fabian (1998) , 
a systematic difference between the electron and the proton temperatures cannot explain the 
inconsistency, and a real departure from hydrostatic equilibrium must happen in some cases. Even 
disregarding these extreme cases, once the temperature profiles are available it is worth: i) to 
check how far from the cluster center the hydrostatic condition can be assumed, specially if strong 
temperature gradients are present; ii) to estimate the correction to the total mass and baryon 
fraction implied by non-isothermality. Moreover some clusters show an anomalously high value 
of the dynamical parameter (3 spec , possibly suggesting deviations from the equilibrium conditions 
and requiring a detailed analysis of the velocity distribution. The galaxy cluster A 2319 has been 
extensively studied in the past, so that many galaxy redshifts are available, ROSAT PSPC images 
can be retrieved from the public archive and a temperature profile based on ASCA data has been 
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published by |Markevitch 1996| (see however [Molendi 1998] |Molendi et al. 1999] ). 

In the present work we combine these data with the F band photometry of galaxies ( Treves^ 
et al. 1992| ), and compare the gas and galaxy density distributions. We generalize the definition of 
the morphological parameter fixo to verify the hydrostatic equilibrium conditions. The analysis 
suggests the consistency of the hydrostatic model in the presence of a temperature gradient. Thus 
we discuss the mass distribution as obtained by adopting a polytropic model or a simple parabolic 
representation of the temperature profile, and we compare the resulting baryon fraction with the 
limits provided by the standard nucleosynthesis calculations, then deriving constraints on the large 
scale baryon segregation and the cosmological parameter J7 Q . 

We use H Q = 50 /150 kms~ 1 Mpc~ 1 . 



2. The galaxy distribution 

The galaxy cluster A 2319 has been studied by several authors in the radio , optical and 



X-ray bands. It is classified as a BM type II-III and as a richness 1, RS-type cD cluster (see Abell 
Corwin and Olowin 1989| ), The galaxy velocity dispersion a ~1800 km s _1 is particularly high. 



However, Faber and Dressier (1977), on the basis of 31 galaxy spectra already suggested that 
A 2319 is actually two clusters nearly superimposed along the line of sight: the main component 
A 2319 A with an average redshift u^=15882 and a velocity dispersion a a = 873l^4g km s™ 1 and 
the second component A 2319 B, located about 8' NW with ^=19074 and <tb = 573^^g km s _1 . 

More recently Oegerle, Hill, and Fitchett (1995) (OHF) measured several new redshifts, applied 
the "<5-test"of Dressier and Shectman (1988) to locate the A and B components, and assigned the 
139 galaxies of known redshift to the component A and B (or to the background/foreground) on 
the basis of their position and redshift, empirically trying to keep gaussian the velocity distribution 
of A 2319 B. They found N,4=100 and N#=28 galaxies in the two components with ^=15727 km 
s , ci4=1324 km s _1 ub=18636 km s" 1 ctb=742 km s _1 respectively. 

To assign individual galaxies to the A and B components we adopted the results of OHF to 

obtained a first order estimate of the cluster positions, average radial velocities ti", and velocity 

(i) 

dispersions cru\ , and we computed the relevant core radii Rc of the two components, where i=A,B 
identifies the component. Then we assumed the following probability distributions of galaxies 
respect to radial velocities v and projected distance b from the relevant cluster center: 
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9i(v) = ■ ^- exp 



where Nu\ the first order estimate of the number of galaxies of the relevant component, and the 

(i) 

bmax is the radius of the circle containing the Ni observed galaxies. This simple parameterization 
is independent of any assumption about the distance along the line of sight of the two clusters 
and their relative motion. Each galaxy is then assigned to the component of higher probability. 
The 99% confidence volumes are also considered for each cluster, and galaxies outside these 
volumes are assigned to the background or foreground. We obtain the new values Nyi=96 and 
N B =24 v a = 15891/cms- 1 , a A = 1235 ± 90kms~\ v B = 18859/sms- 1 , a B = 655 ± 97fcma _1 . 
Since the resulting velocity distributions of the two components do not show strong deviation 
from gaussian distributions, the reported uncertainties a aA and cr aB , on a a and a B respectively, 
have been computed as a 2 . = a 2 ^/[2(N^ — 1)] i=A,B. The effect of membership uncertainty 
can be evaluated as follows. Given a sample of N galaxies with average velocity v and velocity 
dispersion a, the addition of k galaxies with velocity v = v + 5 produces a new velocity 
dispersion a' 2 = (Na 2 + k5 2 )/(N + k — 1). Therefore, to increase a a by more than 2a aA it is 
necessary to include in the sample A more than k=2 galaxies with a recession velocity exceeding 
VB + 2(TB ^va + 3.46<ta . 

Although on the sole basis of the "5 test" there is a 10% probability that A2319B is 
not a physical association, the strong clustering of large 5 values in a region (see OHF fig. 5) 
corresponding to enhanced X-ray emission suggest that it is a physical entity. Moreover, the 



analysis of bound orbits of the two components A 2319 A and A 2319 B led Oegerle, Hill, and 



Fitchett 1995| to the conclusion that "there is a reasonably high probability that these clusters 
are not bound and will never merge". The latter conclusion is supported by the discussion of 
FGB who compare the X-ray images with simulations of cluster collisions (see section 3). The 
above considerations suggest and legitimate the assumption, adopted in the following, that the 
two clusters are separate entities. 

We have added to the spectroscopic information the F band photometry of A 2319, obtained 



by Frevese et al. 1992 from microdensitometric scans of a Palomar 48 inch Schmidt plates, as 



part of a systematic study of the morphology and luminosity functions of galaxy clusters (Trevese, 



Cirimele and Flin 1992| , Flin et al. 1995| , Trevese, Cirimele and Appodia 1996| , Trevese et al. 1997). 



Due to the low galactic latitude (b~ 13°) the field of A 2319 is very crowded and the automatic 
star/galaxy classification is difficult. Thus we have revised the classification and recovered some 
misclassified object. 

To focus our attention on the main component A, we reduced the effect of the B component 
excluding from our sample all the galaxies classified as B (or background/foreground). We 
adopted a fixed center (a = 19 /l 21 m 11.8 s , 5 = +43 56'39"( J2000)), derived from the centroid of 
X-ray emission as computed in a small (2 arcmin) circle around the intensity peak. This point is 
identified with the center of a spherical structure which we assume to represent the A component. 
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We chose a magnitude limit tuf = + 2=16.33 mag and the resulting fraction of galaxies 
without measured redshift is 0.23. Thus, the fraction of galaxies without measured redshift and 
belonging to the B component is of the order 0.23Nb / (N^ + Ng), i.e. 4%, and is not expected 
to affect significantly the galaxy density profile. We fitted with a maximum likelihood algorithm 
the unbinned galaxy distribution using both a King profile a ga i{b) = oo • (1 + (b/r c ) 2 )~ K + o\, 
and with a de Vaucouleurs profile cr ga i(b) = ctq ■ exp(— 7.67(b/r v )' y s'j + cr b5 where the background 
counts Oft, r v , j g , r c and (3 are free parameters, while oo is determined by the normalization to 
the total number of observed galaxies and b is the projected distance from the cluster center. The 
Kolmogorov-Smirnov (KS) test has been applied in both cases and the results are reported in 
Table 1, where the errors reported represent one-sigma uncertainties derived from Monte Carlo 
simulations described in section 4, and Pks(> D) is the probability of the null hypothesis that 
deviations larger than D are produced by random noise. 

The surface density and the fitting profiles are shown in Figure 1. In this case the King 
profile has a slightly higher probability and will be adopted in the following to derive the 
volume distribution by numerical inversion. However, the differences between the two fitting 
curves, specially for b > 0.1h^ Q Mpc, cannot affect significantly the subsequent discussion of the 
hydrostatic equilibrium conditions. 



To obtain the total luminosity of the cluster, Lt ot (r), we fitted with a |Schechter (1976) 



function the unbinned luminosity distribution excluding the brightest member, using a maximum 



likelihood algorithm, adopting a constant a = —1.25 and M* as a free parameter, as in Trevese 
Cirimele and Appodia 1996| . This gives: 



L tot (r) = 10 -o.4(M^ 2 8.43) + r (2) 

r(l + a, ^p) Jo 

where L tot is expressed in units of 10 13 Lq, r in kpc and Ln m is the limiting luminosity 
corresponding to the magnitude limit {Mf = —21.56 mag) of the galaxy sample adopted to derive 
Pgai- The value of M* changes by less than 1% considering a sample which includes the galaxies 
of A 2319 B. The galaxy mass M ga i{r) is then obtained from the total luminosity assuming an 



average mass-to-light ratio M/Lr = 3.32 ± O.^Mq/L^q from van der Marel (1991) , adopting 
F ~ R for bright ellipticals (see [Lugger 1989| ) and computing Fq from the relation V-F=0.40 
(B-V) QKron 198CQ - 

To estimate the virial mass we have evaluated the r.m.s. velocity dispersion a r in four 
concentric rings each containing 1/4 of the galaxies of known redshift of A 2319 A. The four values 
are 1148 km s _1 , 1415 km s _1 , 1327 km s^ 1 , 1193 km s _1 with a r.m.s uncertainty of about 200 
km s _1 . Thus in the following we assume a constant dispersion, derived from the entire A 2319 A 
sample, ay = 1235 ± 90km s _1 

The resulting virial mass My = 3vr6 G cj T 2 /2G is My = 2.89 x lO 15 ^ 1 ^©, namely only 2% 
less than the value given by OHF, since the decrease of of is almost entirely compensated by a 
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slight increase of the projected virial radius be = 2(1/6) 1 ( |5arazin 1988 ), which in our case is 
b G = 1.736/iso 1 Mpc. 



3. The gas distribution 

From the ROSAT public archive we extracted the available Position Sensitive Proportional 
Counter (PSPC) images corresponding to two observations of A 2319, on March 1991 and 
November 1992, which cover a 128x128 arcmin 2 field with pixel size of 15x15 arcsec 2 and an 
effective resolution of about 25" FWHM, in the energy band 0.5-2.0 keV. The exposure maps 
( Snowden et al. 1992j , Plucinsky et al. 1993), providing the effective exposure time of each pixel, 



are also available at the ROSAT public archive for each observation. We divided each image for 
the relevant exposure map and combined the two images with weights proportional to the maxima 
of the exposure map (1514.6 s and 3200.8 s respectively). The resulting image is shown in Figure 
2. 

Feretti, Giovannini & Bohringer (1997) (FGB) discuss the radio structure of A2319, which 
shows a powerful radio halo. They also analyze two substructural features in the X-rays, one 
corresponds to the E-W elongation in the very center of the A component, detected in the image 
obtained with the ROSAT High Resolution Imager, and is interpreted by FGB as an evidence 
of a recent merging process, likely providing the energy for the radio halo . This feature is 
confined within the inner 5 arcmin and does not affect the analysis of the hydrostatic equilibrium, 
particularly in the region where the temperature is not constant, i.e. for r > 5 arcmin. In the outer 
region the cluster structure is rather regular, as for most cD clusters, except for an elongation 
in the direction of the B component. According to FGB, the B component is in a pre-merger 
state and has not yet proceed far enough to disturb the bulk of the gas, as can be argued from a 



comparison with the cluster collision simulation by (Schindler and Miiller 1993) . 



Thus we analyzed A 2319 A as a separate entity, as discussed in section 2. As a first 
approximation we ignored the presence of the B component, we assumed spherical symmetry and 
as in CNT, we derived the volume density p ga s( r ) of the gas by both numerical inversion of the 
projection equation, and by fitting the observed surface brightness with a "/?- model" ( pavalicrc 



and Fusco-Femiano 1976| ) 1(b) = Iq[1 + (b/r c ) 2 } 3 ^ +1 / 2 + J 6) obtaining consistent results. The 



results are also consistent with FGB and with a more recent analysis of Mohr et al. (199"9)| . Then 
we applied the same procedure after the exclusion of the northern half (5 > 43°56'39"( J2000)) of 
the image, to eliminate the effect of the B component. We have evaluated the constant background 
value lb = 8.86 x 10 -4 cts s _1 arcmin -2 in the region b > 3/i^q 1 Mpc. The fitting parameters are 
reported in Table 2, together with the central proton density n . The one-sigma uncertainties have 
been evaluated by Monte Carlo simulations, described in section 4. The central proton density is 
obtained from the relation: 

/ [1 + (b/r c ) 2 ]~ W+1/2 = EM ■ ^ A(i/, T)dv (3) 
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where A(v res t,T) is the rest-frame cooling function corrected for the cosmological dimming 
factor (1 + z) and computed with the code of Mewe et al. (1986) with 0.3 solar abundance 



corresponding to ^ = 1.2 , v m in and v max define the observing band in the rest-frame 
and EM = / °° '(^ L )n 2 (l + ^)~^^dl (see [Barazin 1988|) , and correction for absorption with 
Njj = 8.89 • 1(T 20 cm- 2 QStark et gTHgj ) has been applied to I Q . 



Again our results are consistent with the analysis of FGB. In particular we also find a smaller 
value of the core radius r c of the A component when the northern half is excluded to reduce the 
effect the B component. The observed surface brightness profile and the fitting /?-model are shown 
in Figure 3. 

For consistency with the optical analysis, where we eliminated the galaxies assigned to the B 
component, in the following we adopt the fit obtained after the exclusion of the northern half of the 
cluster image. The corresponding gas density profile p ga s{r) = P gas [^ + (fr/Vc) 2 ] -3 ^ 2 is obtained 
assuming a constant temperature equal to the emission-weighted temperature Tx = 10.0 ± 0.7keV 
QMarkevitch 1996| ). The non-isothermality does not significantly affect the results ( |Markevitch et| 



|al. 1996| ), due to the weak dependence of the emissivity on temperature, in the adopted band 
(see next section). The gas density profile has been computed also by a non parametric numerical 
deprojection as in CNT, obtaining consistent results. 



4. Hydrostatic model, mass distributions and baryon fraction 

As pointed out in CNT, it is possible to check the hydrostatic equilibrium condition in a 
non-parametric way, by a direct comparison of the density profiles of gas and galaxies. Assuming 
spherical symmetry, the equilibrium condition implies ( Bahcall and Lubin 1994| ): 



pm p a 2 din p gas (r)/dlnr + dhiT/dlar 

kT ~ dlnp gal (r)/dlnr + dlna 2 /dlnr + 2A ^ ' 

where oy is the radial galaxy velocity dispersion, A = 1 — (f 1 ) measures the anisotropy of the 
velocity distribution and p is the average molecular weight which we assume equal to 0.58 ([Edge 



and Stuart 199 1|) . For constant ay and T, and A = 0, this implies p gas oc p[ 3 al where f3 spec = iir "^ T 



gal 

is a constant representing the ratio between the energy per unit mass in the galaxies and the 
gas respectively (|Cavaliere and Fusco-Femiano 1976| ). In this case it is possible to define the 
morphological parameter (3xo = d In p gas (r) jd In p ga i to be compared with (3 spec which is obtained 
from the spectroscopic observation of galaxies and the gas temperature derived from X-ray 
spectra. As already discussed in CNT, the very existence of a wide range of densities where (5xo 
is constant supports, in many cases, the validity of the isothermal model. However in the presence 
of a temperature gradient, as in the case of A 2319 A, both [3 sp ec and fixo depend on radius and 
the equilibrium equation reads: 

Pspec = Pxo + Pro, (5) 
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where (3to = dlnT/dln p ga i and we still assume that the galaxy velocity distribution is isotropic 
and a r is constant. To check the validity of equation || we used the temperatures obtained 
by Markevitch 199C . We computed (3 spec , fixo an d /?to a t three projected radii of 4, 10 and 



20 arcmin, corresponding to the boundaries of the four rings whose temperatures are given by 
Markevitch 199C. With the constant velocity dispersion a r = 1235 ± 90km s _1 derived in section 2, 



Pspec ranges from 0.86 ± 0.15 in the first point, to 1.52± 0.47 in the outermost point, i.e. the ratio 
between the energy per unit mass in gas and galaxies is close to unity within a central isothermal 
region of about 0.5 h^Q Mpc, while it decreases in the outer regions. 

Figure 4 shows lnp gas , obtained by numerical deprojection, versus In p ga i- From the slope of 
the curve is possible to derive a value of ftxoif) at each radius. 

The values of Pxo( r ) range from 0.59 ± 0.06 in the first ring, to 0.64 ± 0.06 in the outermost 
ring. The one-sigma uncertainties are evaluated from Monte Carlo simulations described below. 
We also computed /?To( r )> which is close to zero in the cluster center and increases up to 0.58±0.44 
at rwl h^ 1 Mpc. The one-sigma uncertainties are derived from errors on the temperature values 
reported in Figure 2 of Markevitch 1996. Figure 5 shows /3 spec (r) as a function of the RHS of 



equation |, ((3 X o + Pro)- 

We stress that the relation between the two quantities on the x and y axes is not automatically 
implied by their definition, since the galaxy density p ga i{r) only appears in the RHS and is 
observationally independent of the quantities in the LHS. 

The value of j3 spec {r) is slightly larger than (3xo + Pro- This could indicate that a is still 
slightly overestimated due to the presence of the background component A 2319 B. However the 
intrinsic statistical uncertainty does not allow this level of accuracy. Furthermore, even small 
deviation from the spherical symmetry, or from isotropy of the velocity distribution, could have 
comparable effects. 

All we can safely say is that there is an increase of [3 spec versus (flxo + Pro), which is 
consistent with a straight line of unit slope. Thus, within the present uncertainties the data are 
consistent with equation [|, namely with validity of the hydrostatic equilibrium, also in the outer 
region of the cluster where the temperature declines. 

The accuracy of this type of check, as applied to a single cluster, is limited by various factors. 
Although future X-ray data will provide much higher signal-to-noise, the uncertainty on galaxy 
density profile is intrinsically limited by Poisson noise on galaxy counts. Moreover subclustering 
and unknowable deviations from spherical symmetry will always produce an uncertainty on the 
galaxy density deprojection. 

Nevertheless, the systematic application of the method described, to all the cluster with 
measured temperature distribution ( Markevitch et al. 1998; ) will likely provide statistical indication 



on the validity of, or the deviation from, the equilibrium conditions in the outer parts of galaxy 
clusters. 
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Assuming that our results indicate the validity of hydrostatic equilibrium, we can derive the 
distribution of the total binding mass Mt(t) of A2319A: 



, x kT ( din p gas (r) dhiT \ 

Mtot{r) = -J^G { dlnr + dlnVj " (6) 

The slopes of the temperature profiles are crucial in establishing whether the non-isothermality 
causes an increase or decrease of the mass estimate, with respect to the isothermal /3-model. 

In fact, from equation |(| indicating with M^°*(r) the mass derived by an isothermal /?-model 
with temperature T{ so t the fractional change in the mass estimate is: 

M tot (r)-Mi°f(r) = T(r)-T lsot T(r) dWT 

^tot t ( r ) Ti so t Ti so t 

where the two terms on the r.h.s. of the equation can be of the same order. 

It is customary to adopt a polytropic gas distributions as the simplest analytic representation 
of a non-isothermal gas in hydrostatic equilibrium. In this case T oc Pg^si with the polytropic 
index 7 ranging from unity to 5/3 for an isothermal or adiabatic equilibrium respectively. The 
polytropic model implies flCowie, Henriksen and Mushotzky 1987 ): 



P 9 as{r) = p° gas [l + (r/r c ) 2 ]- 5 , T(r) = T [l + (r/r c ) 2 ]" a , (8) 
3/3/2 



1 + 77(7- l)/2' 



a = 5(7 — 1) 



where rj = dlne/dlnT and e is the emissivity of the gas, integrated in the adopted band. We 
can fit the temperature profile of A 2319 provided by |Markevitch 1996 , neglecting the effect of 



projection and adopting 77 ~ —0.2, to determine the polytropic index 7. Notice that, due to 
small value of 77, the result is not significantly different if we simply assume 5 = 3/3/2 and rj = 
in the above expressions. The slope \dT/dr\ of the temperature profile reaches a maximum at 
r/r c = —a + Vet 2 + 1 and progressively decreases in the outer regions. 

As a result, for Ti sot = T a , the quantity A(r) is positive only for r/r c < x y = \J 7 1 /" — 1. This 
limit decreases for increasing (3 and 7, e.g. x 7 ~ 0.78 for /3 = 1, 7 = 5/3 and 77 = —0.2, while 
x-y = V e 2 — 1 ~ 2.53 for j3 = 1/3 and 7 = 1. We stress that A(r) is always negative at large radii, 
and this implies an enhancement of the "baryon catastrophe". The best fit value is 7 = 1.091 and 
the resulting temperature profile is shown in Figure 6. 

The quality of the fit is quite poor, and there is a probability of 91% that the deviations 
are non random. The increment of the mass in the central region, with respect to isothermal 
model, and the enhancement of the baryonic catastrophe obtained with a polytropic model are a 
mere artifact of the particular shape of the polytropic temperature profile, which is not a good 
representation of the data, at least in the case of A 2319 A. 
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More specifically, the data seem to indicate an almost isothermal central region and an 
increasing slope of the temperature profile for increasing radius. 

This behavior is consistent with a change of the polytropic index with radius from an 
isothermal (7 = 1) towards an adiabatic (7 = 5/3) hydrostatic equilibrium (see Sarazin 1988). A 



fit with the law T(r) = T a — ar is shown in Figure 6. In this case the probability is P(> x 2 )=0.99. 

Figure 7 shows Mj*f*(r), corresponding to the isothermal model, and M tot (r) as obtained with 
both the polytropic model and the quadratic interpolation of the temperature profile. On the basis 
of the above discussion we assume the latter as the best representation of the mass distribution. 

In the same figure, galaxy and gas masses, M ga i(r) and M gas (r) are also shown. The latter is 
computed both in the isothermal approximation and using the temperature profile: the result is 
only weakly dependent on the temperature changes. 

Mgas(r) is steeper than M ga i(r) and M tot (r), which have more similar slopes. As a consequence 
the gas mass dominates over the galaxy mass at large radii. These results are consistent with 
previous findings of CNT. 

The statistical uncertainties have been evaluated by Monte Carlo simulations of the entire 
reduction process. For the X-ray data, starting from a "/?- model" corresponding to the fitting 
parameters, we generated 500 random sets representing the photon counts in each radial ring 
with Poisson noise, and we extracted 500 random background values, with a standard deviation 
estimated from intensity fluctuation of the surface brightness outside 3h§Q Mpc where the 
background value has been measured. Then we fitted with a "/?- model" each intensity profile, 
obtaining the statistical distribution of the fitting parameters. Finally we extracted 500 values of 



the temperature in each of the four rings corresponding to the data of Markevitch 1996 , with the 
relevant standard deviations, and we fitted the temperature profile with a parabolic law. Then we 
applied to the simulated data the same algorithms applied to real data for the evaluation of the 
mass profiles. This procedure allows to define a one-sigma confidence interval for the gas mass 
M gas and for the total mass Mtot as a function of radius. A similar procedure was adopted for the 
galaxy distribution. Then we extracted, for each cluster simulation, a random value of M/Lf q 
from a gaussian distribution with a mean value and a standard deviation corresponding obtained 



from van der Mar el (1991). The one-sigma confidence intervals are reported as shaded areas in 
Figure 7. 

We define the luminous mass as the sum of the gas mass and the galaxy mass as deduced 
using an average stellar mass-to-light ratio (see section 2): M\ um = M gas + M ga i. The above 
results imply that the dark matter Mjark = M to t — Mi um has a distribution similar to M ga i(r). 
This provides a constraint on the mechanism of galaxy and cluster formation. 

Since an unknown fraction of the dark matter is baryonic, Mi um /Mt t represents a lower limit 
on the baryon fraction /{,. Figure 8 shows this lower limit as a function of radius, as computed in 
the isothermal approximation and taking into account the temperature gradient by the polytropic 



- 11 - 



model and by the quadratic interpolation. We can compare our estimate of /& with the results of 
Mohr et al. (1999)| who gives the f b values f b = 0.213 ± 0.004 and f b = 0.297 ± 0.027, at lh^ Mpc 



and at 1.91/igQ 1 Mpc respectively, the latter distance representing r$Qo, namely the radius within 
which the mean density is 500 times the critical density p cr u = 3H% /8ttG. Our isothermal values 
fb = 0.180 ± 0.017 and fb = 0.252 ± 0.045 are slightly lower, because were derived by a fit of the 
southern part of the X-ray image, to exclude the effect of the B component. 

At 2 h^Q Mpc fb, as computed with the quadratic interpolation, becomes respectively 66 % 
and 56 % of the isothermal and polytropic values, thus mitigating the "baryon catastrophe". 

Assuming /& ~ 0.2 as typical of galaxy clusters, the residual discrepancy between fij, = fb^o 
and the corresponding value derived from nucleosynthesis calculations Q^ ucl ~ 0.076 ± 0.004/igQ 2 
( Buries and Tytlcr 1998| ), can be reconciled assuming Q Q < Qb Ucl / fb ~ 0.4. 



5. Summary and Conclusions 

We have performed a new analysis of the Abell cluster A 2319, and assigned the individual 
galaxies to the A 2319 A and A 2319 B components by an objective criterion taking into account 
both position and redshift. The resulting velocity dispersion of the A component is slightly 
smaller, but consistent with the previous determination of OHF. 

We have obtained photographic F band photometry of the cluster galaxies, which allows us 
to construct the cluster luminosity function and galaxy density profile. 

We have analyzed archival ROSAT PSPC images separating the A component on the basis 
of the optical information, and we have obtained a gas density profile of the A component. The 
result is consistent with recent studies of Feretti, Giovannini and Bohringer (1997), Mohr et al. 
(1999). 



Since, according to Markevitch 1996] , A 2319 A shows a radial gas temperature decrease, we 



have generalized the method introduced by Cirimele, Nesci &: Trevese (1997) CNT, in order to 
check the validity of the hydrostatic equilibrium in the case of a non isothermal gas. 

We have derived the total mass profile Mt t( r ) through the non isothermal hydrostatic 
equation adding new evidence in favor of the results of CNT that the total mass and the galaxy 
mass have similar radial distributions, more concentrated with respect to the gas component. 

Polytropic models imply smaller masses at large radii, respect to the isothermal model, i.e. a 
higher value of the baryon fraction. Thus, the use of a polytropic model would enhance the baryon 
catastrophe. 

However we have shown that the polytropic model is inconsistent with the observed 
temperature profile, at least in the specific case of A 2319 A. A parabolic representation of 
the temperature profile gives, instead, a total mass larger than computed in the isothermal 
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approximation, mitigating the baryon catastrophe. 

In any case, fb is larger than fl^ ucl ~ 0.076 ± 0.004/igQ 2 , resulting from nucleosynthesis 
calculations and recent measures of the deuterium to hydrogen ratio (D/H), in high resolution 
studies of the Lya forest ( [Buries and Tytler 1998| ). Under the assumption that the value of 
fb ~ 0.2, found for A 2319 A, which is consistent with other estimates ( [Evrard 1997 , CNT, Ettori 



fc Fabian 1999] , Mohr et al., 1999) is typical of galaxy clusters, it is possible to derive the following 



conclusions. If the f2 Q = 1 assumption is kept, then the baryonic fraction within galaxy clusters 
is not representative of the cosmic value and clusters must be surrounded by dark matter halos 
( White and Fabian 1995 ). If, on the other hand, beyond r ~ 2/i^q 1 Mpc the material has not yet 



fallen into the cluster, as infall models indicate, then fb is representative of the cosmic value and 
the cosmological parameter must be f2 Q = Q^ ucl / fb ( [White et al. 1993 and refs. therein), i.e. in 



our case £l Q ^ 0.4, where the uncertainty on the latter value depends on the intrinsic spread of the 
baryon fraction, in the presently observable volumes around galaxy clusters. 

Recently Markevitch et al. (1998) have collected temperature profiles for 30 galaxy clusters 
based on ASCA observations. However, only a minority of clusters is regular enough, specially in 
the outer regions, to allow an X-ray and optical check of the equilibrium conditions as suggested in 
the present paper. This limits the accuracy and the reliability of a "measure" of the cosmological 
parameter Q Q , based on a comparison between the baryon fraction fb and Q^ ud . Moreover future 
cosmic microwave background experiments are expected to provide tighter constraints on f2 D 



( Mandolesi et al. 1995 ) as compared with the results derived from the analysis of the matter 



distribution in galaxy clusters. 

On the basis of this new estimate of Q , once the equilibrium conditions of the non- isothermal 
regions are verified by the systematic application of the analysis outlined in the present paper , 
it will be possible to extend the (otherwise questionable) estimates of the general distribution of 
luminous and dark matter based on the hydrostatic model, to the outer regions of galaxy clusters 
of a statistical sample, thus providing new constraints for the models of cluster formation an the 
physics of large scale baryon segregation. 

We are grateful to the anonymous referee and to the editor for comments and suggestions. 
This work has been partly supported by Ministero dell'Universita e della Ricerca Scientifica e 
Tecnologica (MURST). 
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Table 1. Galaxy density fits. 





de Vaucouleurs profile 
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King profile 
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Table 2. X-ray brightness fits. 



Image 


Jo x 10 4 


n a x 10 3 


r c 


P 


x> 




cts s~ 1 arcmin~ 2 


cm~ 3 


Mpc 






Overall" 


801±144 


5.28±0.70 


0.213±0.090 


0.518±0.070 


174/62 


Southern 3, 


950±171 


6.62±0.84 


0.159±0.071 


0.511±0.069 


82/62 



a For both images exposure time is 4715.4 s and the background is h = 8.86 • 
1CF 4 cts s~ x arani-nT 2 
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Fig. 1. — Radial profile of the surface density of galaxies. Vertical bars represent the standard 
deviation of the number counts in the relevant annulus. The curves represent the fits of unbinned 
data. Solid line: King profile; dotted line: de Vaucouleurs profile. 

Fig. 2. — ROSAT/PSPC image of A 2319. The image is filtered with a gaussian of 100 arcsec 
(FWHM). The contour levels, starting from 0.40 counts pixel -1 (lpixel = 15" x 15") increase by a 
factor 2 1 / 2 . Notice the asymmetry in the N-NW direction, due to the B component. 

Fig. 3. — X-ray surface brightness profile of A 2319 A, after the exclusion of the northern half of 
the image. Errors represent the standard deviation of photon counts. The solid line is the fit with 
the /3-model. 

Fig. 4. — The gas density p ga s( r ), derived from numerical inversion of the projection equation (see 
CNT), versus the fitted galaxy density p ga i(r)- Vertical lines show the relevant radii, indicated in 
Mpc (h$o = 1). The local slope of the curve defines the function (3xo{ r )- 

Fig. 5. — fispec versus (Pxo + Pro) computed at the radii of 4, 10, 20 arcmin, corresponding to the 
boundaries of the four rings whose temperatures are given by Markevitch (1996). 

Fig. 6. — Temperature versus radius of A 2319 A, as deduced from Markevitch (1996). dotted: best 
fit polytropic model; short dash: quadratic interpolation T(r) = T a + ar 2 . 

Fig. 7. — The mass as a function of radius for A 2319 A. Solid: total binding mass M^ t '(r), for 
the isothermal -model with Ti sot = 10 ± 0.7keV; dotted line: M tot (r) for the polytropic model; 
shaded area: M tot (r) for T(r) = T Q — ar 2 ; long dash: galaxy mass M ga i(r); dot dash: gas mass 
M gas (r). Shaded areas represent the la uncertainty as deduced by Monte Carlo simulations. 

Fig. 8. — Baryon fraction /&(r) as a function of radius. Solid: isothermal "/3-model"; dotted: 
polytropic model; shaded area: model with T(r) = T Q — ar 2 with the la uncertainty as deduced by 
Monte Carlo simulations. 
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